Evidence for a Massive Black Hole in the SO Galaxy 
NGC 4342 



Nicolas Cretton 

Sterrewacht Leiden, The Netherlands 

Frank C. van den Bosch 

University of Washington, Seattle, USA 

Abstract. We have constructed axisymmetric dynamical models of the 
edge-on SO galaxy NGC 4342: simple two-integral Jeans models as well as 
fully general, three-integral models using a modified version of Schwarz- 
schild's orbit superposition technique. The two-integral models suggest a 
black hole (BH) of 3 or 6 x 10 8 M©, depending on the data set. The three- 
integral models can fit all ground-based and HST data simultaneously, 
but only when a central BH is included. Models without BH are ruled 
out at better than 99.73% confidence level. We determine a BH mass of 
x 10 8 M Q . This corresponds to 2.6% of the bulge mass, making 
NGC 4342 one of the galaxies with the highest BH mass to bulge mass 
ratio currently known. 



1. The data 

NGC 4342 is a low luminosity (Mb = —17.47) edge-on SO galaxy in the Virgo 
cluster, displaying both an outer disk (r > 5") and a stellar nuclear disk (r < 1"). 

The data consists of long slit spectroscopy obtained from the ground, and 
of multi-color photometry and single aperture spectroscopy obtained with HST 
(see van den Bosch, Jaffe & van der Marel 1998). We modelled this complex 
mass distribution using the Multi Gaussian Expansion technique of Emsellem, 
Monnet & Bacon (1994). 

2. Jeans models 

We have first constructed simple dynamical models assuming a distribution func- 
tion (DF) of the form DF = f(E,L z ). Under the assumption of constant mass 
to light ratio (M/L), we have solved the Jeans equations and compared the pro- 
jected velocity dispersions (after seeing convolution and pixel binning) to the 
observations. Although this model shows some discrepancies with the data, it 
suggests a BH mass in the range (3 - 6)xlO 8 M (see Cretton & van den Bosch 
1998 for details). These models cannot fit simultaneously the two data sets, 
suggesting a more sophisticated type of modeling (three- integral DF). 
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Figure 1. Fits to the kinematic data for NGC 4342: mean streaming 
velocity V, velocity dispersion a and Gauss-Hermite moments I13 and 
/14 for ground-based long slit data (William Herschel Telescope) and 
HST FOS apertures. Model A (dotted line) has no BH, model B (solid 
line) a 3 x 10 8 M BH and model C (dashed line) a 6 x 10 8 M BH. 

3. Three-integral models 

Our scheme follows Schwarzschild's method with an extension towards the use 
of kinematic constraints: We first deproject the surface brightness to obtain a 
three-dimensional mass density distribution (assuming a constant M/L). From 
the Poisson equation the gravitational potential corresponding to the stellar 
body is derived. At that stage, a dark component (BH or dark halo) can be added 
to the potential. In this total potential, a large library of orbits is computed and 
projected onto the observable space (position on the sky, line-of-sight velocities). 
Again the necessary steps of seeing convolution and pixel binning are performed 
for a fair comparison to the data. Finally the orbital weights are computed that 
best reproduce the photometric and kinematic constraints (see Figure 1 for the 
kinematic fits of 3 models). 

The inclusion of the complete distribution of line-of-sight velocities greatly 
helps in constraining the final model. The three-integral models have more free- 
dom than the simpler Jeans models and have no problems fitting small features 
of the data (e.g. in the rotation curve). Moreover, contrary to the two-integral 
models, they can fit simultaneously both data sets (ground-based and HST). 

To determine the exact values of M/L and BH mass, we compute a grid 
of models in (M/L, BH) and use x 2 ~statistics to assign confidence levels (see 
Figure 2). On the basis of the complete data set, models without BH can be 
ruled out at a confidence level > 99.73 %. The best fit values for a (M/L, BH) 
are {G-3±oi M /L , 3.0t\i x 1O 8 M ): model B in Figure 2. Therefore in this 
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Figure 2. Contour plots of iso-x 2 levels. The first three contours de- 
fine the formal 68.3, 95.4 and 99.73 (thick contour) per cent confidence 
levels. Solid dots indicate actual model calculations. The % 2 surface 
of panel a results from excluding h$ and /14 as well as all HST/FOS 
measurements from the constraints. Panel b shows the x 2 plot when all 
constraints (HST and ground-based) are used. The asterisks labeled 
A, B and C indicate special models discussed in Section 4. 

galaxy, the ratio BH mass to bulge mass (MBH/-^buigc) is the highest currently 
known (2.6%), together with the case of NGC 3115 (see e.g., Ho 1998). 

4. Internal kinematics 

Figure 3 shows the second order velocity moments and velocity dispersions of 
models A, B and C averaged over a cone of opening angle 15° around the equa- 
torial plane. This region contains the inner and outer discs, but in the very 
center, the bulge contribution dominates the light. 

Models A, B and C mainly differ in the center (due to different BH masses). 
Outside 3", (w 2 ) dominates the dynamics in accordance with the rapid rotation 
of the outer disk. From 3" to 12", we pass from the azimuthally anisotropic 
bulge into a radially anisotropic region: the dynamically cold outer disk built 
up of close-to-circular orbits with low a^. 

Models B and C have a z /&R remarkably constant at ~ 0.9, and are thus 
not too different from two-integral models (for which this ratio is exactly 1.0). 
The same has been observed in M32, one of the few galaxies for which three- 
integral models have been constructed (van der Marel et al. 1998). Merritt 
and collaborators (see this volume) have proposed that, under the influence of 
a central BH, box orbits are destroyed and the global mass distribution evolves 
towards axisymmetry (see also Gerhard &; Binney 1985). At the same time, this 
tends to erase any dependence on a third integral in the DF, consistent with what 
we have found here. When MBH/-^buige is > 2.5%, the very short (~ one crossing 
time) evolution towards axisymmetry provides a negative feedback mechanism 
that limits the BH mass by cutting off its fuel suply (Merritt & Quinlan 1998). 
Therefore in this scenario, the maximum BH mass that can be accreted in the 
center is 2.5% Mb u i gc , close to the observed maximum. 

The main change going from model A to model C, is a strong increase 
of (J<p/(JR in the inner ~ 3". In this region the circular velocities increase 
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Figure 3. Internal kinematics in the equatorial plane for models A, 
B and C. The first lines shows the second moment and the bottom line 
the normalized velocity dispersion cr a /a tota ,\- The insert in the upper 
right panel explains the line convention. 

strongly with increasing BH mass. Nevertheless, all three models provide an 
almost equally good fit to the observed rotation velocities because, going from 
model A to C, the ratio of — L 2 -orbits over +L 2 -orbits increases in the interval 
0.5" < R < 3.0". This causes the net streaming motions of all three models to 
be roughly similar despite the large differences in circular velocities and explains 
the strong increase in u^. 
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